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Abstract 

The so-called Vegetation Red-Edge (VRE), a sharp increase in the reflectance 
around 700 nm, is a characteristic of vegetation spectra, and can therefore be used as a 
biomarker if it can be detected in an unresolved extrasolar Earth-like planet integrated 
reflectance spectrum. Here we investigate the potential for detection of vegetation spec- 
tra during the last Quaternary climatic extrema, the Last Glacial Maximum (LGM) 
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and the Holocenc optimum, for which past climatic simulations have been made. By 
testing the VRE detectability during these extrema when Earth's climate and biomes 
maps were different from today, we are able to test the vegetation detectability on 
a terrestrial planet different from our modern Earth. Data from the Biome3.5 model 
have been associated to visible GOME spectra for each biome and cloud cover to derive 
Earth's integrated spectra for given Earth phases and observer positions. The VRE is 
then measured. Results show that the vegetation remains detectable during the last 
climatic extrema. Compared to current Earth, the Holocene optimum with a greener 
Sahara slightly increases the mean VRE on one hand, while on the other hand, the 
large ice cap over the northern Hemisphere during the LGM decreases vegetation de- 
tectability. We finally discuss the detectability of the VRE in the context of recently 
proposed space missions. 

Keywords: Extrasolar planets, Earth-like planets, vegetation, biomarker, biosigna- 
ture, Earth spectrum, modern Earth, Holocene optimum, Last Glacial Maximum. 

1 Introduction 

One of the most challenging goals of the future space missions is to detect life on extrasolar 
planets. These missions are designed to image a planet as a single dot floating beside its 
parent star, using high-contrast imaging techniques. Photons from the unresolved planet 
will feed a low-resolution spectrograph and give first insights into planet's chemistry. Re- 
cent works suggest that Earth's vegetation is detectable in the Earth integrated spectrum 
at visible wavelengths [Arnold et al. 20021 IWoolf et al. 20021 lArnold 2008] , assuming ex- 
traterrestrial plants, if they exist, are similar to Earth's plants. Our green vegetation 
indeed has a reflectance spectrum showing a sharp increase - the so-called vegetation red 
edge (VRE) - around 700 nm (Fig. [I]) that could be considered as a possible global 
biomarker. 
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A question arising from these results is whether vegetation was also detectable when 
the (early) Earth looked significantly different than today. In other words, if an Earth- 
like planet with a climate different from the modern Earth's climate is detected one day, 
will we be able to detect vegetation on it? To answer this, we investigate the vegetation 
spectral detectability during the last Quaternary climatic extrema, for which past climate 
simulations have been made by general circulation models (GCMs). This simulated climate 
is used as a driver for an equilibrium vegetation model, providing biome maps of the Earth 
at these epochs. 

The Last Glacial Maximum (LGM) occurred about 21 kyr (kilo-year) ago. Tempera- 
tures globally of the order of 4°C colder than today [Braconnot et al. 2007] were responsible 
for a large extent of sea and continental ice sheets, especially in the northern Hemisphere. 
Sea level was also significantly lower (—121 ± 5 m [Fairbanks 1989]) resulting in larger 
emerged lands. 

During the Holocene optimum (6 kyr before present, hereafter kyrBP), Earth's north- 
ern Hemisphere was ~ 0.5°C warmer than today [Braconnot et al. 2007] . The sea level 
was rising, although still slightly lower than today, because the amount of water released 
from continental ice was partially compensated by an isostatic response of the continents, 
i.e. a rise of the continents following the deglaciation of the northern Hemisphere. The 
Sahara was more vegetated than today |Ritchie &: Haynes~1 987, Bonfil s~et al. 2001] . its 
desertification being the main large-scale change in land cover for the last 6000 years 
[Glaussen et al. 1999j . 

Here, based on a combination of biomes and simulated cloud cover maps for the peri- 
ods mentioned above, associated with 300 to 800 nm spectra for each biome and cloud, we 
build Earth integrated spectra, for given Earth phases and observer positions. The vege- 
tation spectral signature is then extracted from these spectra. The first part of the paper 
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Figure 1: Reflectance spectra of photosynthetic (green) vegetation, non-photosynthetic 
(dry) vegetation and soil (from [Clark 1999J). The so-called vegetation red edge (VRE) 
is the green vegetation reflectance strong variation from ~ 5% at 670 ran to ~ 70% at 
800 ran. Note that the spectrum of a non-vegetated soil also shows a higher reflectance in 
the near-IR than in the visible. 

describes the model of the Earth from which we derive global reflectance spectra. The 
second part of the paper discusses the results. Note that some results showing the rotating 
Earth and corresponding instantaneous spectrum are available on-line as animated gif files 
([http://www.obs-hp.fr/~arnold/results/2009JJAstrobio/VRE . html]) . 



2 Main inputs for the Earth model 

2.1 Biomes maps from Biome3.5 model for 0, 6 and 21 kyrBP 

The equilibrium terrestrial biosphere model BIOME3 [Haxeltine Prentice 19 96] is a pro- 
cess based vegetation model simulating biogeography and biogeochemistry. Model inputs 
consist of monthly climate variables (temperature, precipitation, and sunshine) data, min- 
imum annual temperature and atmospheric C02 concentration. A coupled hydrological 
model calculates moisture availability on the basis of a soil texture class. The model simu- 
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lates the productivity and cover (Leaf Area Index) of a set of plant functional types. The 
comparative dominance of these types is then used to estimate the dominant biome for 
a given grid cell. A total of 24 biomes may be simulated, including barren ground and 
snow/ice cover (Fig. [2]and[3j). 

The model was run for the three periods of interest: modern, Holocene optimum 
and LGM Earth (Fig. [2] and [3]). For the modern period, gridded climate data on a 
grid was taken from [Leemans Cramer 199T] , and atmospheric C02 concentra- 
tion was set to a pre-industrial level of 280 ppm. For the two past periods, climate data 
was obtained from simulations run for under the Paleoclimate Modelling Intercompari- 
son Project- We used data from the UK Meteorological Office General Circulation Model 
(UKMO GCM), run under climatic forcings for these two periods [Hewitt fc Mitchell 1 996 , 
IHewitt fc Mitchell 1997] . The GCM output was downscaled to the same resolution grid as 
the modern period by calculating anomalies for each past period. The GCM control run 
value was subtracted from the past value for each GCM grid point, and then these anoma- 
lies were added to the modern climatology grid points. Atmospheric C02 concentration 
was set to 265 ppm and 180 ppm for the Holocene optimum and LGM respectively. 

The results obtained agree well with observational evidence of past vegetation change 
[Prentice et al. 2"0 00], The mid-Holocene is dominated by a greening of the Sahara, and 
a poleward shift of the high latitude treelines. The LGM is characterised by a large 
expansion of steppe and grassland vegetation at mid-latitudes, with a near disappearance 
of temperate forest. There is also the fragmentation of tropical forest and widespread 
expansion of tundra vegetation at high latitudes. There are, however, some disagreements 
with the data, for example the expansion of grassland in Australia. 

1 Paleoclimate Modelling Intercomparison Project (PMIP), |http:/ /pm ip. lsce.ipsl.fr/ 
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Figure 2: Map legend and biomes index. 



2.2 Biome spectra 



We used a collection of GOME spectra taken the 18, 19 and 20th of May 1999, corre- 
sponding to springtime and the middle of the greenup period in the northern Hemisphere. 
In May the average NDVI (Normalized Differential Vegetation Index) is just half of its 
maximum value reached in Jul}o Sun declination is ~ +19.7° for the observing dates in 
May. It is fall in the southern Hemisphere and the polar night is present at latitudes below 
—70.3°. The GOME spectra are reflectances obtained from the ratio of a given ground 
pixel radiance R over solar irradiance I©, 



it R 



Iq cos 6 s 



(1) 



2 http:/ /earthobservatory.nasa.gov/Newsroom/NewImages/Images/PIA04335.gif 
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Figure 3: Maps of biomes today (top), 6 kyrBP during the Holocene optimum (middle) and 
21 kyrBP during the LGM (bottom), with Continental and sea ices. Note that 6 kyrBP, 
a large part of the Sahara displayed a temperate grassland biome and almost no barren 
regions, while today and at 21 kyrBP, it consists of barren and desert (shrubland and 
steppe) landscapes. 



where 9 S is the solar zenit angle. A reference spectra for each biome r& is extracted from 
the GOME data (Fig. 0]). In order to find a given biome spectrum with minimal cloud 
contamination, the algorithm choses the spectrum with the lowest signal in the interval 
450 to 500 nm, considering that the signal in this spectral range rapidly increases with the 
presence of clouds. This algorithm fails for three types of surfaces: ocean, snow/ice and 
cloud. These spectra are chosen 'by hand'. 

It must also be noted that this algorithm chooses the most vegetated pixel for a given 
biome: if a cloudless pixel is a mixture of sand and vegetation, it will be ignored in 
favor of a pixel with more vegetation, which will display a lowest signal in the spectral 
interval mentioned. This suggests that the reference pixel (and spectrum) for each biome 
is probably biased toward the most vegetated pixel of a given biome area. 

Lastly, since it is spring in the northern Hemisphere and fall in the southern Hemisphere, 
the reference spectrum for the biomes of deciduous species will be preferentially picked up 
in the 'greener' northern Hemisphere. Therefore, our simulations will be more relevant 
for reconstruction of the northern Hemisphere rather of the southern Hemisphere, with 
a probable overestimate of the vegetation signal if the southern Hemisphere is mainly in 
view. 

2.3 Cloud cover 

2.3.1 Modern cloud cover reflectance from ISCCP data 

To build a realistic Earth's spectrum, we need both the cloud cover map and the reflectance 
map of that cloud cover, r cc . We used cloud covers from the International Satellite Cloud 
Climatology Project (ISCCPg). We selected only data for the month of May from the 1999 
to 2004 Vis-IR observations, to be consistent with the GOME spectra. Since low, middle 



3 http:/ /isccp. giss.nasa.gov/ 
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Figure 4: Biomes spectra measured by GOME the 18, 19 and 20th of May 1999. The 
spectra have been recorded at air-mass (here total path length through the atmosphere) 
ranging from 2 to 2.6 depending on solar and spacecraft line-of-sight angles. See Fig. [2] for 
the biome legend. Note that all plots do not have the same vertical scale. 



Table 1: Quantitative data about cloud cover normalized reflectance from real data (IS- 
CCP) and model (biome3.5). 





Modern Earth 


Modern Earth 


Holocene optimum 


LGM 




(ISCCP) 


(biome3.5) 


(biome3.5, 6 kyrBP) 


(biome3.5, 21 kyrBP) 


mean R cc 


0.46 


0.40 


0.40 


0.39 


max R cc 


0.86 


1.00 


1.00 


1.00 


min R cc 


0.00 


0.03 


0.03 


0.05 


+60 ' +90 ° 


0.51 


0.50 


0.50 


0.41 


+40 ° +60 ° 


0.53 


0.52 


0.52 


0.58 


-40 ° +40 ° 


0.38 


0.22 


0.22 


0.24 


-60 ° -40 ° 


0.62 


0.81 


0.81 


0.75 



and high altitude clouds have different typical optical depth and reflectances, we build the 



cloud cover reflectance r cc from the sum of the three different altitudes cloud cover maps 



ccl, ccm and cch, weighted by their albedo, at,, au and an respectively 0.69, 0.48 and 
0.21 [Mana be &; Stric kler 1964], and normalized to the albedo of low altitude clouds cll, 



i.e. 



rcc = cc L H cc M H cch- (2) 

CLL CLL 



The result from this formula is that the three different types of clouds contribute to 



a given pixel's reflectance proportionally to their cover over that pixel. The formula also 



can lead to a few pixels with reflectance above 1, which are then simply clipped to 1. 



The obtained map (Fig. [5]) is normalized to the albedo of low clouds for which we have a 



reference spectrum from the GOME data. 



2.3.2 Cloud covers from Biome3.5 model for 0, 6 and 21 kyrBP 



Estimates of past changes in cloud cover were taken from the UKMO GCM, as used as 



input for the Biome 3.5 model. The and 6 kyrBP cloud covers are the same (Fig. |SJ). The 



simulated cloud cover at 21 kyrBP shows locally overestimated cloud cover in the southern 
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Figure 5: Normalized reflectance map of the cloud cover estimated from ISCCP data on 
low, mid and high cloud cover maps and clouds albedo jManabe &: Strickler 1964] for the 
mean of cloud cover for month of May for six years from 1999 to 2004. No data are available 
for the Antarctica because this region is always in the night in May. Quantitative data are 
given in Table [TJ 
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Figure 6: Normalized reflectance maps of the cloud cover estimated from cloud cover 
simulated by the UKMO GCM, with the same grey-scale as for Fig. [5l for modern Earth 
or at 6 kyrBP (top) and 21 kyrBP (bottom). In the 21 kyrBP map, artifacts induced by 
the existence of a wavetrain resulting from the passage of the lower atmospheric circulation 
over the topography of South America are visible in the southern Hemisphere: The cloud 
cover is slightly overestimated in the middle of the Atlantic Ocean and at the South of the 
Madagascar island. Quantitative data are given in Table [TJ 
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Hemisphere, caused by the existence of a wavetrain resulting from the passage of the lower 
atmospheric circulation over the topography of South America (Fig. [6]). The simulated 
cloud cover output consists of 17 layers of clouds at given pressures from 1000 (ground) 
to 10 mbar (highest clouds). The model shows that the highest layers, at 30, 20 and 10 
mbar, are very high clouds above Antarctica which are therefore not relevant for our study. 
The remaining 14 layers are quadratically summed into three layers of low, middle and 
high clouds, according the ISCCP classification versus cloud pressurqj (i.e. low clouds for 
pressures from 1000 to 680 mbar, middle clouds from 680 to 440 mbar, and high clouds 
above 440 mbar). A simple sum of the 14 layers into 3 would mean that the clouds in each 
layer contribute linearly and cumulatively to the total cloud cover over a given pixel. This 
is not realistic and leads to overestimation of cloud cover (many pixels have a cloud cover 
significantly larger than 100%). Using a quadratic sum allows the super imposition of layers 
to be taken into account, leading to much more realistic results when the kyrBP model 
is compared with real ISCCP data. The cloud reflectance maps are calculated following 
Eq. [2] as done with the ISCCP data. 



3 Disk-averaged spectra from the model and extraction of 
the vegetation signal 

A 2-layer map of the Earth is created with biome and cloud cover for each pixel. For 

ocean, 2-layers are used too, but the biome is replaced with the sea ice cover. Once this 

composite map is built, we proceed with an orthographic projection where the central 

longitude and latitude corresponds to the position chosen for the virtual distant observer. 

A realistic image of the Earth -i.e. a radiance map-, with limb darkening and phase, is 

finally obtained by multiplying the previous projection by a mask of a white lambertian 
4 http:/ /climate. msrc.sunysb.edu/cpt/isccp-cltypes.gif 
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sphere at the desired phase. We consider indeed a lambertian Earth, i.e. all biomes, ocean 
and clouds have an isotropic reflectance. Thus no specular reflection on ocean surface is 
taken into account. Specular reflection has indeed an albedo of 0.02, moreover attenuated 
by lands and clouds [Breon &: Henriot 20 06]. We estimate from [Breon Henriot~2 006j 
that specular reflection would increase the overall Earth albedo by ~ 0.01 and therefore 
has a negligible quantitative impact. 

We assume that a lambertian approximation of biomes reflection, with respect to real 
anisotropic BRDF, has negligeable impact on the global spectrum. Specular reflection of 
the Sun in ocean would make the Earth brighter and consequently decrease the relative 
signal of vegetation. Our model also considers all clouds as optically thick, according to 
ISCCP mean optical thickness of clouds^ equal to 3.86. We understand that to derive more 
accurate spectra, our model would probably benefit from radiative transfer calculations, 
an implementation of the sun glint effect and the BRDFs of all scattering surfaces included 
in the model (biomes, oceans, clouds). 

The global spectrum is computed from the projection of the Earth described just above: 
The global Earth Reflectance spectrum ER(\) is a linear combination of the spectra of 
each biome, ocean (ice-free or partly covered by sea ice) and clouds, weighted by their 
respective number of pixels and integrated radiance in the projected image, a pixel being 
itself a combination of two or three elements (biome and cloud, or sea, sea ice and cloud). 
The result is normalized to the flux integrated in the image of a white lambertian sphere 
observed with the same phase, giving the result spectrum in terms of apparent albedo A v 
|Qiu et al. 2003 . 

We want to measure the height of the VRE embedded in Earth reflectance ER(A), 
considered as an example of a large-scale ground-based signature of life. We want to 
5 ISCCP data analysis http://isccp.giss.nasa.gov/climanall.html 
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remove all a priori known biases, these biases being due to atmospheric molecules of biotic 
or abiotic origin, or due known minerals spread over large land surfaces. The VRE is 
quantitatively approximated by the ratio [Arnold et al. 20021 IHamdani et al. 2006J 

VRE = ri ~ rR (3) 

where 77 and rjj are the near-infrared (NIR) and red reflectances integrated over spectral 
domains defined on both sides of the Red Edge. This formula gives correct results only 
if the spectrum has been cleaned from known sources of atmospheric bias, as described 
above. There are at least three sources of bias that can be corrected a priori and that have 
been considered here. 

The first is the absorption of atmosphere molecules. Ozone has a weak but wide ab- 
sorption over the full visible range (Chappuis band), and absorption is non-zero at the 
wavelengths of interest for the VRE. The Earth spectrum is thus divided by a reference 
ozone absorption spectrum, extracted from a GOME spectrum of white clouds and adjusted 
to the data. Water vapor and oxygen absorptions are also corrected with MODTRA^jfl 
reference spectra: the band depth of the reference spectra are optimized to minimize the 
residual of the fit, after which the spectrum is divided by the fit. These two molecules 
have lower impact on the VRE extraction than ozone because their absorption bands are 
relatively narrow, moreover the spectral domains over which the VRE is measured can be 
chosen between H2O and O2 absorption bands. 

Scattering in the atmosphere is the second source of bias we take into account. Rayleigh 
and aerosol scattering produces a negative slope in the spectrum in the blue part of the 
spectrum. We consider here that the aerosol and the Rayleigh contribution to scattering 
are equal at A e = 700 nm [Lena 1 996], 

6 MODTRAN is distributed by Ontar Corporation. 
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The third source of bias we take into account is the spectral signature of deserts that 
adds a smooth positive slope across the visible and near-IR domains where the VRE is mea- 
sured [Arnold et al. 2003} lArnold 2008] . and the function defined to fit Rayleigh, aerosol 
and desert biases from 350 to 685 nm is thus written 

= + A4+ ^.5 A 1.5 + A * >< >< *of ■ (4) 

The spectrum is our model of desert reflectance, extracted from a GOME spectrum 
taken over Sahara and corrected for Rayleigh and ozone biases, and so 3 is the spectrum 
of ozone extracted from a GOME spectrum of white clouds. The ER(A) spectrum is 
then normalized to the vegetationless model /(A) to measure the VRE using the spectral 
domains [650:685nn] and [750:758nm]. It can be noted here that the correction of H2O 
and O2 is only relevant for the fit function between 580 and 680 nm. 

To test the robustness of these bias corrections, we built an Earth model where all 
pixels of continents were desert, with oceans, clouds, sea and continental ices remaining 
unchanged. We obtained a mean VRE of 0.1% with a standard deviation of 0.2% after an 
integration over a full Earth rotation, for an observer seeing mostly northern Hemisphere 
of an half-phase Earth (inclinasion i = 60 ° , orbital phase <f> = 180 ° - eastern elongation). 
For identical geometrical parameters, the vegetated Earth model gives a VRE of 5.9% 
(standard deviation of 2.9%) suggesting that the /(A) function well fits a vegetationless 
spectrum, i.e. meaning that the biases mentioned above are well corrected to allow the 
detection and the measurement of VRE values as low as « 1%. 
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Table 2: 24h-averaged disk-averaged VRE (%) and apparent albedo A v measured from 
our Earth model during springtime for different Earth phases and observer's positions. 
The angle <p is the orbital phase: = 0° is for Western maximum elongation, = 90 ° 
for superior conjunction, <j> = 180 ° for Eastern elongation and <f> = 270 ° for inferior 
conjunction. 
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4 Results and discussion 



The reconstructed scenes of the Earth, from which the results below are extracted, have 
been made using a Sun declination of 5 = +19.7°, corresponding to the Earth-Sun con- 
figuration for May. The Earth phase viewed by the observer is reconstructed for several 
inclination i. The angle <p is the orbital phase of the Earth viewed from the location of the 
observer: = 0° is for Western maximum elongation, <fi = 90 ° for superior conjunction, 
<fi = 180 ° for Eastern elongation and 4> = 270 ° for inferior conjunction. 

In Table [5J the VRE and albedo values obtained from the model are given for the 
different periods. For the modern Earth, the albedo values obtained with the ISCCP 
data are higher than with Biome3.5, due to a lower cloud cover with Biome3.5, leading 
to a lower VRE. This suggests that the Biome3.5 model probably overestimates the VRE. 
Nevertheless, relative comparisons between Biome3.5 results remain consistent although 
the absolute values may be overestimated. 

The main result from this study is that the climate differences between the Quaternary 
extrema and the modern Earth's climate have little impacts on the vegetation signal, since 
the spectral signature of vegetation is not washed out during the LGM and does not increase 
very much during the Holocene optimum. There are nevertheless interesting differences in 
'vegetation visibility' between these two extrema and the situation today, depending of 
the observer's position: When the northern Hemisphere is mainly in view, vegetation is 
indeed more visible during the Holocene optimum as shown in Fig. [7] and Table [2j The 
VRE was typically 6%, thus higher than today (5%) or during the LGM (4%). This is 
mainly due to a greener Sahara as shown in Fig. Eta). This is true for an observer with 
mostly the northern Hemisphere in view, i.e. inclination % = 45 ° . In inferior conjunction 
configuration, the crescent-shaped Earth mainly shows its northern part. There are no 
significant differences between modern and Holocene Earths (VRE ~ 2%) and vegetation 
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becomes almost undetectable for the Earth with the LGM large polar ice cap (< 1%). If 
the observer is above the northern pole (observing from ecliptic North pole) with in view 
a half-illuminated Earth, there are again no significant differences between modern and 
Holocene Earths (VRE ~ 4%). However, as the North pole ice is much more extended 
during the LGM, the signal of the vegetation is significantly weaker (only 2%), due to the 
higher reflectance of a larger polar zone and the consequently smaller vegetated areas is 
view. 

The three bottom lines from Table [5] give results for Earth views in which southern 
Hemisphere contribution is high, with inclination i ranging from 90 to 180 ° . As pointed 
out in Section 12.21 our model probably overestimates the VRE values of scenes with the 
southern Hemisphere mostly in view as our calibration was made during fall time (May) 
in that hemisphere and deciduous species should have a lower VRE than with the spectra 
used in the simulation that have been collected in the northern Hemisphere. Nevertheless, 
there are two trends to be mentioned: the first one is that the difference between Holocene 
and modern period becomes negligible, because the Sahara contribution decreases or even 
vanishes from these southern view points. The second remark is that, paradoxically, the 
southern Hemisphere VRE during the LGM is higher: this is induced by a greener Australia 
during the LGM, as we obtained a simulated wider distribution of temperate grassland and 
less desert in Australia using the Biome3.5 model (Fig. [3|) at the Holocene optimum. As 
discussed above, this does not agree with observations of past vegetation and may arise 
from an incorrectly simulated climate over this region. 

Figures [10] and [TT] finally show modern Earth's scenes with corresponding spectrum 
obtained following the method described in Sect. El In Fig. [T0| we show the spectrum for 
central longitudes of —20 ° (Africa and Europe) and 120 ° (Northern Pacific Ocean and 
Alaska) at which the model gives the maximum and minimum VRE, respectively of 8.1 
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and 1.5%. Fig. [TOl show the Earth spectra based on Biome3.5 data for central longitudes 
of —20 ° (Africa and Europe) at which the model gives the maximum VRE, i.e. 8.6, 11.0 
and 6.3% at 0, 6 and 21 kyrBP, respectively. 

Will next space observatories be able to detect VRE-like features in an Earth-like planet 
spectrum ? The Terrestrial Planet Finder Coronagraph mission (TPF-C) is aimed to make 
the direct detection and the characterization (multi-spectral imaging at resolution up to 
70) of nearby terrestrial extrasolar planets in the visible range. The telescope would have 
an elliptical aperture of (6 to 8) x 3.5 m 2 |Stapelfeldt et al. 2005[ IShaklan et al. 2006] . An 
Earth-like planet in the habitable zone is 10 10 fainter than its parent star, making such 
planets very faint objects in the V=28 to 31 range. According to |Stapelfeldt et al.~2 005 



exposure time of the order of one day will be needed to detect them and spectral character- 
ization will add 'weeks of additional integration time per target'. More quantitatively, in 



a comparison of different types of coronagraphs |Guyon et al. 2006 , the authors conclude 
that an Earth-like planet at about 10 pc could be detected in 2.3 h with a signal-to-noise 
ratio of SNR=7 with a 100 nm bandwidth from 0.5 to 0.6 \im with a 4-m telescope, and 
even in 14 minutes with a 6-m telescope, both equipped with one of best performers of 
the comparison, a phase-induced amplitude apodisation coronagraph. Now if we want to 
be able to detect a 4% VRE-like feature in a spectrum at a resolution of 70 around 700 
nm, the SNR should be about 80 for a 3-sigma detection, or 125 for a 5-sigma detection. 
With the 4-m coronagraphic telescope metioned above, the total exposure time becomes 
vey long, 18 and 44 weeks for 3 and 5-sigma detection, respectively. This in agreement with 
the order of magnitude given by [Stapclfcldt et al. 20 05 . With a 6-m telscope, the total 
exposure times decrease to 13 days and one month, respectively. If a larger telescope than 
the TPF-C-like telescopes mentioned above is considered, the exposure time will decrease 
again. According to |Guyon et al. 2 006 , a planet could be detected in 2 minutes with 
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a 12-m telescope. The total exposure times to detect a 4% feature in a spectrum would 
then decrease to 2 and 5 days, for 3 and 5-sigma detection, respectively. These very long 
exposure times suggest that features of a few percents will be extremely difficult to detect, 
maybe even not within reach of a 'first-generation' 4-m TPF. A 6-m telescope seems to be 
the minimum to envisage the detection of seasonal variations. With exposure times of 2 to 
5 days, a 12-m telescope would probably not freeze the rotation of the planet and rather 
give an integrated spectrum of the planet over several rotations, but atmospheric changes 
at times scales of the order of a week or longer would be possible to detect and follow up, 
as well as longer seasonal variations. 

Unfortunately TPF-C might not fly before 2025, and in the meantime, a smaller in- 
strument could look for bigger and closer terrestrial planets: A 2-m class space corona- 
graph like SEE-COAST could indeed detect a VRE on super-Earths closer than 3 to 6 
pc [Schneider et al. 2008] . depending on their radius whose maximum value is presently 
unclear. These super-Earths are as well habitable than Earth-sized planets. 

More ambitious instruments that might fly later this century will allow us to make 
a step furthermore in the characterization of these planets and the observation of 'exo- 
vegetation'. Hypertelescopes in space - interferometric sparse arrays of small telescopes - 
will indeed allow us to see Earth-like planets as small resolved disks several resels across 
|Labeyrie 1999 Labeyrie et al. 2008 and this clearly will help us to detect photosynthetic 
life on these planets! For example, a 150-km hypertelescope would provide 40 resolution 
elements (resels) across an Earth at 3 pc in yellow light |Labeyrie~l 999|. And a formation of 
150 3-m mirrors would collect enough photons in 30-min to freeze the rotation of the planet 
and produce an image with at least ~ 300 resels, and up to thousands depending on array 
geometry (Fig. 1 12|) . At this level of spatial resolution, it will be possible to identify clouds, 
oceans and continents, either barren or perhaps (hopefully) conquered by vegetation. 
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Figure 7: Vegetation Red Edge measured with Biome3.5 data for clouds, ice and biomes, 
and GOME spectra (and ISCCP for clouds at kyrBP only). The VRE calculated with 
the ISCCP clouds and the BIOME3.5 model are in correct agreement, although the latter 
is typically 1% higher. Central longitude is decreasing, according to the eastward Earth 
rotation. 




Figure 8: Vegetation Red Edge measured with Biome3.5 for clouds, ice and biomes, and 
GOME spectra (and ISCCP for clouds at kyrBP only). The VRE calculated with the 
ISCCP clouds are significantly lower than those calculated from the BIOME3.5 model, 
probably due to the difference between the corresponding cloud maps for the southern 
Hemisphere as suggested by Table [Q Central longitude is decreasing, according to the 
eastward Earth rotation. 
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Figure 9: Earth scenes reconstructed for the observer's positions given in Table at 0° of 

central longitude. Scenes are for A = 750 nm, a wavelength at which vegetation is brighter 

than in the visible. For a better readability, the radiance is shown here at power 0.7 to 

enhance the darker regions of the scene and underline the shape of the sunlit Earth seen 

94 

by the observer. 



Figure 10: Earth scenes and corresponding spectra reconstructed for two observer's posi- 
tions at inclination i = 45 ° (Tableland Fig. [7|), for central longitudes of —20° (top, 
Africa and Europe) and 120 ° (bottom, Northern Pacific Ocean and Alaska) at which the 
model gives the maximum and minimum VRE, respectively. The Earth spectrum is in 
white, and a vegetationless model in yellow (moreover with O2 and H2O bands corrected) 
allows to visualize and calculate the VRE. For the record, cloud and ocean spectra are 
shown in blue. Scenes are for A = 750 nm, a wavelength at which vegetation is brighter 
than in the visible. For a better readability, the radiance is shown here at power 0.7 to 
enhance the darker regions of the scene and underline the shape of the sunlit Earth seen 
by the observer. 
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Figure 11: Earth scenes and corresponding spectra reconstructed for one observer's position 
at inclination i = 45 ° (Tableland Fig. [7]), for central longitudes of —20° (Africa and 
Europe) at which the model gives the maximum VRE. The Earth spectrum is in white, 
and a vegetationless model in yellow (moreover with O2 and H2O bands corrected) allows 
to visualize and calculate the VRE. For the^ecord, cloud and ocean spectra are shown in 
blue. Scenes are for A = 750 nm, a wavelength at which vegetation is brighter than in 
the visible. For a better readability, the radiance is shown here at power 0.7 to enhance 
the darker regions of the scene and underline the shape of the sunlit Earth seen by the 
observer. 



Figure 12: Simulated image of the Earth at 3 pc (10 light-years) observed with a 150-km hy- 
pertelescope interferometric array made of 150 3-m mirrors working at visible wavelengths 
(from |Labeyrie 1999| ). North and South America are visible. Note that this simulation is 
done at visible wavelengths, while in the (very) near-IR at 750 nm, vegetated areas would 
be much brighter and more easily detectable on continents. Spatial resolution at 750 nm 
would remain the same than at visible wavelength with the same hypertelescope flotilla 
being spread over 225-km instead of 150-km. 

5 Conclusion 

In order to mimic an Earth-like extrasolar planet observed without spatial resolution, i.e. 
as a single dot, in a different climate state than our contemporary Earth, we considered the 
Earth during the last climatic extrema, 6 and 21 kyr ago, and studied the vegetation spec- 
tral signature integrated over the full illuminated phase seen by the observer. The results 
presented in this paper suggest that vegetation disk-averaged signature remained visible in 
Earth spectra during these climatic extrema, with VRE values close to those for the mod- 
ern Earth. Our discussion emphasized some differences induced for example by the greener 
Sahara during the Holocene optimum, or by the large northern ice cap during the LGM. 
Results would probably benefit in accuracy from radiative transfer calculations, an imple- 
mentation of the sun glint effect and the BRDFs of all scattering surfaces included in the 
model. The spectra would thence also display a more accurate Rayleigh scattering 'slope' 
and absorption bands depth. Further, improved GCM estimates for these past periods are 
now available, including coupled ocean-atmosphere simulations [Braconnot et al. 2007] . as 
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well as dynamic vegetation models |Sitch "e£ al. 2003] . both of which may help improve the 
past vegetation distributions obtained here. All these improvements could be considered 
for future developments. 

Of course a next step in this work could also be the consideration of the Earth at much 
older epochs, when continents had a very different shape than today and climates maybe 
much more extreme than those considered in this work, although paleoclimatic models 
are probably less confident for these epochs. It would although probably provide us with 
relevant illustrations of Earth spectrum variability over one full day, and related VRE 
variations. Current climate change, thence biome maps for the end of this century, might 
also tell us if vegetation will remain a detectable biomarker for observers out there. 
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